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ABSTRACT 

Because of the semi-coUisional nature of the solar wind, the coUisionless or exospheric approach 
as weh as the hydrodynamic one are both inaccurate. However, the advantage of simphcity makes 
them useful for enlightening some basic mechanisms of solar wind acceleration. Previous exo- 
spheric models have been able to reproduce winds that were already nearly supersonic at the 
exobase, the altitude above which there are no collisions. In order to allow transonic solutions, a 
lower exobase has to be considered, in which case the protons are experiencing a non-monotonic 
potential energy profile. This is done in the present work. In this model, the electron velocity 
distribution in the corona is assumed non-thermal. Parametric results are presented and show 
that the high acceleration obtained does not depend on the details of the non-thermal distribu- 
tions. This acceleration seems, therefore, to be a robust result produced by the presence of a 
sufhcient number of suprathermal electrons. A method for improving the exospheric description 
is also given, which consists in mapping particle orbits in terms of their invariants of motion. 

Subject headings: acceleration of particles — methods: numerical — stars: winds, outflows — solar wind 



— Sun: corona 

1. Introduction 

Most cosmic bodies eject matter into space, 
but the solar wind is the first and, up to now, 
only stellar outflow to have been measured in situ 
(Neugebauer 1997). Numerous sophisticated mod- 
els have been developed since Parker's pionner- 
ing papers (Parker 1958, 1960) using complicated 
mechanisms. The different acceleration mecha- 
nisms, the origin of the high-speed solar wind as 
well as the associated problem of coronal heating 
have been recently reviewed in a comprehensive 
way by Cranmer (2002). However the solar wind 
acceleration and its properties are still far from 
being well understood. The reason of this diffi- 
culty is that the solar wind is neither a collision- 
dominated medium nor a collisionless one. The 
Knudsen number K„, which is defined as the ratio 
of the particle mean free path and the density scale 
height, is close to unity at Earth's orbit (see e.g. 
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Hundhausen (1972)). This means that neither the 
hydrodynamic approach nor the pure collisionless 
one (also called exospheric) are fully appropriate 
to model the solar wind expansion and to explain 
its observed properties. 

The classical fluid approach is applicable for the 
extreme regime when Kn <^ 1, that is when the 
medium is collision-dominated. In this case, the 
particle velocity distribution functions (VDFs) are 
Maxwellians as the medium is assumed to be at 
local thermodynamic equilibrium. The Euler or 
Navier-Stokes approximations are applicable and 
produce a thermally driven wind out of the hot 
solar corona. There are two problems with this 
approach. Firstly, the particle VDFs might not 
be Maxwellians at the base of the solar wind. 
Secondly, the particle VDFs are not Maxwellian 
in the solar wind. There is an increasing num- 
ber of both theoretical (Vihas et al. 2000; Leub- 
ner 2002), and observational evidences (Esser and 
Edgar 2000) , which tend to show that non-thermal 
VDFs can develop and exist in the high corona 
and even in the transition region. This is because 
in a plasma the particle free paths increase rapidly 
with speed (cx w*), so that high energy tails can de- 
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velop for Knudsen numbers as low as 10"'^ (Shoub 
1983), i.e. even in collisional plasmas. A fortiori, 
high energy tails can be expected to be found in 
the weakly collisional corona and solar wind ac- 
celeration region. Indeed, it is well known that 
the solar wind electron VDFs permanently exhibit 
non-thermal tails that can be modelled by a halo 
Maxwellian population (for e.g. Feldman et al. 
(1975)) or by the power law part of a generalized 
Lorcntzian or Kappa function (Maksimovic ct al. 
1997b). In the frame of the fluid approach, which 
intrinsically cannot handle suprathermal tails, the 
effect of non-thermal VDFs on the solar wind ac- 
celeration can be understood through an increase 
of the heat flux (Hollweg 1978; Olbert 1981). 

An alternative way of taking into account the 
possible effects of coronal non-thermal distribu- 
tions is to use a kinetic approach. Among the 
various kinetic approaches for the solar wind, the; 
simplest one is probably the exospheric one, which 
totally neglects binary collisions between particles 
above a given altitude called the exobase. The 
first solar wind model of this type was developed 
by Chamberlain (1960) by analogy with the evap- 
oration of planetary atmospheres. This first ex- 
ospheric model, modelling the radial expansion 
of the solar corona from the thermal evaporation 
of hot coronal protons out of the solar gravita- 
tional field, produced a solar breeze. The subsonic 
speed obtained by the theory was partially due to 
an inadequate assumption: the electrostatic field 
was taken so as to ensure hydrostatic equilibrium 
(Pannekoek 1922; Rosseland 1924), which is in- 
consistent with an expanding atmosphere. The 
improved exospheric models by Jensen (1963) and 
Brandt and Cassinelli (1966) were the first to be 
able to reproduce supersonic solar wind flows. In 
these models, multiple exobase locations were as- 
sumed, in order to take into account the vari- 
able mean free path of the particles as a func- 
tion of their velocity. However these models still 
used the inadequate Pannekoek- Rosseland electric 
field as an imposed external solution. The actual 
(outward) ambipolar electric field, which ensures 
plasma quasi-neutrality and zero electric current 
is greater, thereby accelerating protons to greater 
speeds. Models using this correction (Lemaire and 
Scherer 1971a; Jockers 1970) produced supersonic 
winds, but too small speeds for explaining the fast 
solar wind (~ 700 - 800 km s'^). 



More recently, Maksimovic et al. (1997a) have 
generalized these calculations by considering non- 
Maxwellian velocity distribution functions for 
the electrons in the corona, e.g. generalized 
Lorentzian or Kappa functions. With such non- 
Maxwellian distributions having suprathermal 
electron tails, a higher electrostatic potential is 
needed to ensure zero charge and current, there- 
fore producing larger terminal bulk speeds. In 
essence, this comes about because the electron 
tail tends to increase the escaping electron flux, so 
that, to preserve quasi-neutrality, the electrostatic 
potential increases in order to trap more electrons, 
which in turn accelerates the protons outwards. 
This model yields a reasonable description of bulk 
solar wind properties, giving densities, tempera- 
tures and speeds within the ranges observed at 
1 AU, even though the details of the VDFs are 
not reproduced, as expected since collisions are 
neglected. Its major interest is the prediction 
of high speeds without assuming extremely large 
coronal temperatures and/or additional heating 
of the outer corona, as is needed in hydrodynamic 
models. More basically, the main achievement of 
exospheric models is to fiirnish a possible driving 
mechanism for the fast solar wind, with a single 
assumption: the suprathermal electron VDF at 
the exobase. 

However, the Maksimovic et al. (1997a) model 
cannot be applied for low altitude exobases, which 
is the case of coronal holes, from where emanates 
the fast solar wind. In these deep coronal lay- 
ers, the gravitational force acting on the protons 
is stronger than the electric one, so that the total 
potential for the protons is attractive out to some 
distance where the two forces balance each other. 
Farther out, the outward electric force dominates. 
This means that the total potential energy for the 
protons is not monotonic, presenting a maximum 
at a certain distance from the exobase (Jockers 
1970) and therefore not all the protons present 
at the exobase are able to escape. The presence 
of such a maximum is not taken into account in 
the Maksimovic et al. (1997a) model, nor in the 
Lemaire and Scherer (1971a) one, since both mod- 
els started beyond the vicinity of the sonic point. 

The purpose of the present paper is to study the 
efifect of non-thermal electron VDFs in the frame 
of a transonic exospheric solar wind model. We 
use a special technique described by Jockers (1970) 
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and Khazanov ct al. (1998) consisting in mapping 
particle orbits in terms of the invariants of motion. 
The same problem has been recently considered 
by Lamy ct al. (2003) using a different formula- 
tion, which involves an approximation relative to 
the escaping particles rate (see Appendix A). The 
present work, which uses the same approximation, 
sets the basis of an exact exospheric description of 
the solar wind acceleration. The reader should, 
however, have in mind that the present model is 
still a very simplified one, which does not pretend 
to describe all solar wind properties. Instead, it 
may be useful to determine some basic aspects and 
a possible driving mechanism of the solar wind, 
avoiding ad hoc assumptions on energy dissipation 
and using as few free parameters as possible. 

In section 2 we recall the basics of an exo- 
spheric solar wind model. In section 3 we outline 
the difficulties arising when dealing with a non- 
monotonic potential energy for the protons and 
describe the technique used to calculate the in- 
terplanetary electrostatic potential. In section 4 
we deal with non-Maxwellian electron distribution 
functions and consider three cases: a Kappa, a 
sum of two Maxwellians (a core and a halo) and 
a sum of a Maxwellian and a Kappa. The latter 
distribution is rather general and reproduces the 
main features of VDFs observed in space plasmas: 
a Maxwellian profile at low speeds and a high en- 
ergy tail with a power law shape. In section 5 
we describe the results of our model and compare 
them with observations. A summary and final re- 
marks are given in section 6. 

2. Basics of a kinetic coUisionless model 

In exospheric or kinetic coUisionless models, a 
specific altitude, called the exobase, is defined as 
an abrupt boundary between the collision domi- 
nated region (in which a hydrodynamic approach 
is valid) and a completely coUisionless one. This 
boundary is usually defined as the distance rg from 
the center of the Sun where the Knudsen num- 
ber Kn is equal to unity, i.e. where the Coulomb 
mean free path becomes equal to the local den- 
sity scale height. Although the mean free path 
depends on the particle temperature and, there- 
fore the exobase should be different for protons 
and electrons, a common exobase is assumed in 
order to simplify the calculations. Furthermore, 



in order to concentrate on the basic physics and 
since the exobase location is defined from param- 
eters that are not accurately known, we put the 
exobase at the surface of the Sun, i.e. ro = IRq. 
This stands as an approximation for a more real- 
istic value that should be slightly greater and has 
no important qualitative impact on the results. 

The coUisionless nature of the plasma above 
the exobase allows one to calculate the velocity 
distribution of each particle species as a function 
of the distribution at the exobase by using Liou- 
ville's theorem with conservation of energy E and 
magnetic moment /x in order to solve the Vlasov 
equation for a time stationary wind. These two 
constants of motion are defined as: 

o 

Tnv 

E = + m(j)g + Ze(j)E (1) 



where v is the speed of the particle of mass m, Ze 
its charge, (j)g{r) = —MqG/v the gravitational po- 
tential, 4>Eir) the interplanetary electrostatic po- 
tential and v± the velocity component perpendic- 
ular to the magnetic field vector B. The choice of 
using the constants of motion {E, ii) instead of the 
speed and pitch angle is convenient as it removes 
the spatial dependence in the distribution func- 
tion. In a purely coUisionless model, the magnetic 
moment is conserved, even though this may not be 
true for protons in the high-speed wind (Schwartz 
and Marsch 1983). 

The classification of particles in different species 
is based on their trajectories along magnetic field 
lines. These trajectories depend on the velocity 
and pitch angle at the exobase, or correspondingly 
on their energy and magnetic moment. There are 
four different classes of trajectories: inc;oniing, es- 
caping, ballistic and trapped particles, as defined 
in Lemaire and Scherer (1971a). Escaping are the 
particles that have enough energy to escape from 
the Sun, while ballistic are those with insufficient 
energy that are returning towards it. Another kind 
of non-escaping particles are the trapped ones, 
which do not have enough energy to escape, but 
whose inclination to the magnetic field lines is 
large enough that they are reflected by the mag- 
netic mirror force before reaching the exobase ro- 
Finally there are particles coming from infinity, 
called incoming, which are neglected. Indeed, due 
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to the postulated absence of collisions above the 
exobase, no particles, in principle, are backscat- 
tered in the downward loss cone. This assumption 
can be relaxed in future applications of our model, 
but it is not unreasonable in the fully coUisionless 
case. 

Once a VDF /q is assumed at the exobase level 
ro, the VDF / at any larger distance r is uniquely 
determined by Liouville's theorem. The electron 
and proton distributions and fp and their mo- 
ments - in particular the electron and proton den- 
sities He and Up - are functions of the electric 
potential (pEir). The quasi- neutrality condition 
ne(r) = np{r) is then used to determine, by an 
iterative method, the value of the potential at 
any altitude r. The iterative process is stopped 
when the estimation of ^£;(r) is adequate within 
the required precision. In addition we impose the 
equality of the electron and proton fluxes at the 
exobase, in order to ensure a zero electric cur- 
rent. (Note that the current will also be zero 
everywhere because fluxes of both kinds of parti- 
cles vary as r~^). Then, it is possible to calculate 
the other moments, e.g. pressures and heat flux, 
at any distance r. Strictly speaking, one should 
solve Poisson's equation instead of imposing quasi- 
neutrality. However, for scales much greater than 
the Debye length, which is the case here, the quasi- 
neutrality condition provides a very good approx- 
imation for the electrostatic potential distribution 
in the solar wind. 

Another assumption is the use of a radial in- 
terplanetary magnetic field B{r), varying as r^^, 
which means that the rotation of the Sun as well as 
a possible super-radial expansion of the wind are 
neglected. As was shown by Chen et al. (1972) 
with a fluid approach and more recently by Pier- 
rard et al. (2001) with an cxosphcric approach us- 
ing Kappa distributions, a spiral magnetic field 
does not change significantly the wind density and 
bulk speed. However, temperatures are modified 
and especially their anisotropics. In any case the 
magnetic field is not very far from radial up to 
about lAU so that this approximation is reason- 
able in the framework of the present model in 
which we study fundamental aspects of the fast 
solar wind coming from coronal holes, i.e. rather 
high heliolatitudes. Nevertheless, a super-radial 
expansion of the wind is an ingredient that might 
influence the results, but it is neglected in this zero 



order approximation in order to avoid additional 
free parameters. 

3. Non-monotonic proton potential en- 
ergy 

Solar wind protons are subjected to the at- 
tractive gravitational potential 0g(r) of the Sun 
and to the repulsive interplanetary electrostatic 
field 0_E(r). Their total potential energy is thus 
^{r) = m(j)g{r) + e(j)E{r), which reaches a maxi- 
mum at some radial distance Vmax- Such a maxi- 
mum occurs because the electric force is expected 
to decrease slower than r~^, so that this force, 
which pushes the protons outwards, should dom- 
inate gravity at large distances (Jockers 1970; 
Meyer- Vernet et al. 2003). Below this altitude 
fmax, the gravitational force is larger than the 
electrostatic one and the total potential is attrac- 
tive. The opposite is true above rmax, forcing all 
the protons present at these altitudes to escape. 

Until now, exospheric models using suprather- 
mal distributions for the electrons have consid- 
ered a monotonic potential profile for the pro- 
tons. This means that the exobase was located 
above the maximum, i.e. at ro > Vmax^ so that 
the wind was already supersonic at the exobase 
level. All protons were escaping, as they were ex- 
periencing a monotonic repulsive total potential. 
Even though these models gave supersonic solu- 
tions, they were incomplete because they did not 
consider the transition from subsonic to supersonic 
speeds. Note that in the general case the Parker's 
critical point does not coincide with the potential 
maximum r^ax and is generally below it (Meyer- 
Vernet et al. 2003), even though these two points 
are relatively close to each other. This is also true 
in recent kinetic simulations taking into account 
Coulombian collisions (Landi and Pantellini 2003). 
In the special case of an increasing temperature 
with T oc r"^/**, these two points coincide (Scudder 
1996). 

For coronal holes, from where emanates the fast 
solar wind, the relation ro > Tmax is no longer 
valid. Indeed, in that case, the exobase is located 
deeper in the corona, because of the lower density 
of coronal holes. At these altitudes, some protons 
cannot escape from the gravitational well of the 
Sun and become ballistic or trapped. Such a case 
has been considered by Jockers (1970), who, how- 
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ever did not consider non-thermal distributions 
for the electrons. Recently, Liemohn and Khaz- 
anov (1998) have given a more general theoreti- 
cal framework for dealing with arbitrary potential 
profiles, while Khazanov et al. (1998) have applied 
it to describe the precipitation of a magnetospheri- 
cally trapped hot population and the outflow from 
the high-latitude ionosphere. 

With a non-monotonic potential, the validity 
of previous cxospheric models is questionable be- 
cause of the violation of constraints given by Chiu 
and Schulz (1978) regarding the first and second 
derivatives of the electrostatic potential with re- 
spect to the magnetic field. The technique de- 
scribed by Liemohn and Khazanov (1998) consists 
in mapping particle orbits in terms of the invari- 
ants of motion. This technique has the interest of 
removing the spatial dependence of the distribu- 
tions, which can now be written as f{E,ij). The 
density, flux and pressures are then given by the 
following expressions: 



V2nB 



P|l = ttB j y/E-nB- ^f{E, n)dEdii 

(5) 

V2nB^ f lJif{E,,i) 
J s/E- 



3/2 



f{E,lA 



y/E- y^B-^ 



dEdn (3) 
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j f{E,ii)dEdf, (4) 



P± 



zdEdii (6) 



The spatial dependence is now transferred into the 
region of integration (over which the moments are 
calculated) in E — ^ space. The basic problem 
is the accessibility of the different particle pop- 
ulations in this space. This technique has been 
recently outlined by Zouganelis et al. (2003); for 
a full detailed description of the method, we refer 
the reader to the paper by Liemohn and Khazanov 
(1998). In Appendix A, we describe some elemen- 
tary concepts of this technique, as applied to the 
solar wind. 

4. Non-thermal electron distributions 

Maxwellian distributions can be used for the 
solar wind protons, but solar wind electrons per- 
manently exhibit significant suprathermal tails in 



their distribution functions. The fundamental rea- 
son for the electrons to be non-Maxwellian is that 
fast electrons collide much less frequently than 
slow ones because of their greater free path, so 
that they cannot relax to a Maxwellian. Usually 
the solar wind electron distributions are fitted to 
a sum of two Maxwellians: a core of cold electrons 
and a halo of hot ones (see for instance Feldman 
et al. (1975), Pilipp et al. (1987)). An alternative 
to the core/halo model has been proposed by Va- 
syliunas (1968), that is a generalized Lorentzian 
or a Kappa function, which has been introduced 
by Scudder (1992) to explain high coronal temper- 
atures by velocity filtration. An interesting com- 
bination of these two distributions is the sum of a 
cold Maxwellian and a hot Kappa as a halo com- 
ponent. The sum of two Maxwellians is then a 
particular case of this last one. 

4.1. Kappa function 

The electron distributions in the solar wind are 
observed to have important high velocity tails and 
a convenient way to fit them is to use Kappa (or 
generalized Lorentzian) functions. This was done 
by Maksimovic et al. (1997b) using Ulysses data. 
Some recent observations (Lin et al. 1997) suggest 
the existence of a superhalo-Maxwellian popula- 
tion that does not seem to be fitted by a Kappa 
function. However, this part of the distribution 
{v ^ Vth) does not contribute in our model. The 
main interest of Kappa distributions is that they 
require one fewer parameter than the core/halo 
model, whereas having a power-law suprathermal 
tail as often observed in space (Vasyliunas 1968). 
Recently, there have been some attempts to find 
a physical explanation of these distributions (Col- 
fier 1993; Ma and Summers 1999; Treumann 1999; 
Leubner 2002) and to study how non-thermal elec- 
tron tails can be generated in the lower chromo- 
sphere (Vinas et al. 2000) and the corona (Vocks 
and Mann 2003). 

Maksimovic et al. (1997a) have considered for 
the first time Kappa distributions in a kinetic col- 
lisionless model of the solar wind, but with the 
restriction of a monotonic potential energy for the 
protons, thus resulting in a wind that is already 
supersonic at the exobase. The Kappa function is 
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defined as: 



f (.-t\ "eO r(K+l) I , V 



2 \ 



(7rKt;2)3/2 r(K - 1/2) V 

(7) 

where r{x) is the Gamma function and v^, is the 
thermal speed defined by: 



2k — 3 k^Tf, 



K TOe 



1/2 



(8) 



where fc^ is the Boltzmann constant and rUf, the 
electron mass. 

For speeds v smaller or comparable to Vk,^ the 
Kappa distribution, for any value k > 2.5, is 
rather close to a MaxwcUian having the same ther- 
mal speed. However, the equivalent Kappa tem- 
perature T^, (defined from the second moment of 
the VDF, as the ratio between pressure and den- 
sity) is related to the Maxwellian one Tm, by 
= {k/{k - 3/2))Tm. For v > v^, the Kappa 
distribution decreases with w as a power law (/„ a 
^-2(k-i-i)^_ In the limit oo, /^(f) reduces to a 
Maxwellian distribution with T^^oo = Tm- Note 
that when electron distributions measured in the 
solar wind are fitted with Kappa functions, the 
parameter Kg for the electrons ranges from 2 to 5 
(Maksimovic ct al. 1997b). 

The basic equations for an exospheric model 
with Kappa distributions are given by Pierrard 
and Lcmairc (1996). Those concerning the elec- 
trons are still valid in our model since the poten- 
tial energy of electrons is an increasing function 
of the distance. When considering a Kappa VDF, 
we therefore use the equations (5)- (34) by Pierrard 
and Lemaire (1996), which give the density, flux, 
pressure and energy flux of the electrons above the 
exobase; these quantities correspond respectively 
to the zero, flrst, second and third order moments 
of the Kappa VDF. 

4.2. A sum of two Maxwellians 

Kappa VDF is the most convenient function 
modelling distributions with high energy tails. 
However, there are several reasons for consider- 
ing other kinds of non-thermal VDFs. Firstly, 
there is currently no agreed-upon physical expla- 
nation for Kappa VDFs and, even though they 
are rather close to a Maxwellian at low velocities, 
they may not be as close to it as observed. Sec- 



ondly, moments of order higher than 2k — 1 di- 
verge since oc v"^**"^ at large speeds. Note 
however that this is merely a mathematical dif- 
ficulty, which has no impact on our model since 
it does not involve these higher order moments. 
Thirdly, from Liouville's theorem, the value of k, 
which represents the non-thermal character, re- 
mains constant with altitude in the absence of 
collisions. Actually, one should expect that as dis- 
tance increases and there arc less collisions, the 
electron VDFs should present stronger suprather- 
mal tails, so that k should decrease with distance. 
Indeed, recent observations in the corona suggest 
distributions having a non-thermal character that 
increases with altitude (Esser and Edgar 2000). 

For these reasons we also consider a sum of two 
Maxwellians, which is the classical way of repre- 
senting electron VDFs in the solar wind (Feld- 
man et al. 1975; Pilipp et al. 1987). Such a VDF 
does not have any of the above-mentioned dis- 
advantages. If Uco and Tgo are respectively the 
electron density and temperature at the exobase 
for the cold (core) component and Uho and Tho 
are the same quantities for the hot (halo) compo- 
nent, we can define their relative importance by 
two parameters: ao = nho/rico and tq = Tho/Tco. 
The total density and temperature at the exobase 
is then given by: no = nco{l + ao) and Teo = 
((1 + aoTo)/(l + ao))Tco- At any altitude the den- 
sity is just the sum of the two densities, the core 
and the halo one. The same is true for the fiuxes 
and pressures. What is interesting is that now the 
parameters a and r are functions of the distance, 
so that the non-thermal nature of the distribu- 
tion is not held constant contrary to the case of a 
Kappa distribution. Note that in this purely col- 
lisionless model the core and halo components do 
not interact with one another (as in e.g. the hydro- 
dynamic core-halo model of Chen et al. (2003)). 
The density, fiux and pressures are defined as the 
moments of the distribution function for each par- 
ticle population (escaping, ballistic and trapped) 
as given in Pierrard and Lemaire (1996). Ana- 
lytical expressions of these quantities can be cal- 
culated in function of the electrostatic potential 
(pEir). For a Maxwellian VDF, they are similar to 
those given in Lemaire and Scherer (1971b) and 
will not be expressed here. These expressions are 
to be used for both the core and the halo compo- 
nents. 
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4.3. A sum of a Maxwellian and a Kappa 
function 

A more general form for the distribution func- 
tion is the sum of a Maxwehian core and a Kappa 
halo. This function is closer to a Maxwellian at 
low speeds than a Kappa distribution and has the 
advantage that the non-thermal character is not 
held constant as was the case with a Kappa VDF. 
We use the same definitions of ao and tq as be- 
fore. Note that for k — > oo the VDF reduces to a 
sum of two Maxwcllians. The expressions of the 
moments used are those for a Maxwellian for the 
core and the equations (5)- (34) by Pierrard and 
Lemaire (1996) for the Kappa halo. 

5. Results and Discussion 

As explained in section 2, the exobase location 
is approximated to be at tq = IRq with no serious 
impact on the results. We assume for the tempera- 
tures at the exobase Teo = 10^ isT and Tpo = 2Teo, 
in the range of values observed in coronal holes 
(Cranmer 2002). The density at the exobase does 
not affect the results of velocity or temperature 
and is just a multiplicative factor in the density 
profiles. 

Let us first consider a Kappa VDF for the elec- 
trons as described in section 4.1. The calculated 
electric potential and the total potential energy 
of the protons are plotted in Figures 1 and 2 re- 
spectively for different values of k ranging from 
K = 6 to K = 2.5, a case with a conspicuous 
suprathermal tail. Note that we use k > 2 in 
order for the energy flux to be finite. One sees 
that the value of the maximum of potential in- 
creases and its distance v^ax 

decreases as k de- 
creases. This is because with more suprathermal 
electrons, a stronger electric potential is needed to 
preserve quasi-neutrality. For a Maxwellian VDF 
{k — > oo), the total proton potential energy in- 
creases monotonically and tends asymptotically to 
zero (remaining always negative). 

The bulk speed - the ratio between flux and 
density - is shown in Figure 3 for a Kappa distribu- 
tion. A high terminal bulk speed (> 700 km s~^) 
is obtained when the suprathermal tail is conspic- 
uous (k = 2.5). This is due to the large value of 
the maximum in ion potential energy (r^ 14A;bTpo), 
which is transformed into kinetic energy of the es- 
caping protons as they are accelerated above Vmax- 



An important remark is that the major part of 
this high terminal bulk speed is obtained within a 
small heliocentric distance (~ lOi?©); this is due 
to the large acceleration represented by the large 
slope of the potential above rmax- Note that this 
is the largest terminal bulk speed obtained by this 
model with a Kappa VDF. 

The density profiles are shown in Figure 4. One 
sees that they are nearly independent of the value 
of K. The density at lAU depends on the one 
taken at the exobase. In this figure we assumed 
an exobase density no = 1.8 • 10^^ m~^, which 
corresponds to the density of a coronal hole ex- 
trapolated to ro = li?0 as given by Koutchmy 
(1977) in line with recent studies on atmospheric 
and coronal electron densities (Esser and Sasselov 
1999). With this density, the model yields a den- 
sity of about 6 cm~^ at lAU, of the same order as 
observed in situ. Note that the rest of the results 
do not depend on the density. 

This analysis bears out previous results of an 
exospheric model using Kappa VDFs (Maksimovic 
et al. 1997a). There are however two basic dif- 
ferences. In the present work the velocity pro- 
files span the whole domain from the subsonic 
to the supersonic regime, which was not the case 
when the exobase was located above rmax- The 
second difference is that we obtain high bulk 
speeds with more reasonable temperatures at the 
exobase. Note however that a direct compari- 
son cannot be done because of a slightly differ- 
ent exobase definition resulting in different proton 
temperatures at the exobase. In any case both 
models can produce high bulk speeds without as- 
suming an additional (ad hoc) heating mechanism 
in the outer corona, as is generally postulated in 
hydrodynamic models. Furthermore, the fact that 
a faster wind is obtained with low values of n 
agrees with observations showing that VDFs have 
large suprathermal tails in the fast solar wind but 
are closer to a Maxwellian in the slow wind (Mak- 
simovic et al. 1997b). 

The large suprathermal tails for low values of k 
have another important consequence. They make 
the electron temperatures increase considerably 
with distance up to a maximum (w 7x10^ K) 
within a few solar radii. This maximum in elec- 
tron temperature is smaller for larger values of k 
and disappears as k ^ oo {Fig.5) as does the 
maximum in the total potential energy of protons. 
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This temperature increase is a direct consequence 
of filtration of the non-MaxweUian VDF by the 
attracting electrostatic potential (Scudder 1992). 
This large temperature increase is not observed, 
which suggests that Kappa functions may not be 
adequate to model VDF having suprathermal tails 
in the corona. 

Let us now consider the results obtained with 
electron distributions made of a sum of two 
MaxweUians or a sum of a Maxwellian core and a 
Kappa halo. On the whole the results are rather 
similar. For the same acceleration we obtain ap- 
proximately the same temperature increase as 
in the Kappa case as we can see in Figure 6 
{ao = 0.03, To = 5 and k = 2.5), which corre- 
sponds to a terminal bulk speed of ^ 770 km s~^. 
We deduce that the temperature increase is not an 
artefact of Kappa VDFs, but a general behavior 
of non-thermal distributions. For a given termi- 
nal bulk speed the filtration mechanism results in 
the same temperature increase. It is important 
to note that coUisionless models are expected to 
give correct electron temperatures (Meyer- Vernet 
and Issautier 1998), because collisions with other 
particles do not significantly affect the electron 
energy, whereas collisions between electrons do 
not change their total temperature. In any case 
we should remind that the present model is still 
a zero-order one, which is intended to explore the 
basic physics of the wind acceleration, but should 
not be expected to reproduce all observations in 
a detailed way, because it involves very few free 
parameters. 

Figure 6 shows also the contributions of the dif- 
ferent particle species to the total electron tem- 
perature. At large distances the temperature pro- 
file is the sum of a term oc r^'^/^ plus a constant. 
The r^^/^ term comes from the isotropically dis- 
tributed electrons (ballistic and trapped) confined 
by the heliospheric electric potential, which is 
found to have the same radial variation at large 
distances. The constant term comes from the par- 
allel temperature of the escaping electrons. This 
agrees with analytical results by Meyer- Vernet and 
Issautier (1998) that do not depend on the VDFs 
in the corona, but were obtained with a monotonic 
proton potential profile. When the proton poten- 
tial energy is non-monotonic, the asymptotic elec- 
tron temperature profile is still the sum of a term 
varying as r~^^^ plus a constant (Meyer- Vernet 



et al. 2003); but the relative importance of these 
terms is not necessarily the same. 

In Figure 7 we show results for a sum of two 
MaxweUians. The diagram shows contours of the 
terminal bulk speed (at lAU) as function of ao 
and To, where we can see that this kind of VDF is 
able to explain the values of the fast wind speed 
( 700 — 800 km s~^). The terminal bulk speed in- 
creases with increasing To, which is not surprising 
since the halo temperature increases. Concerning 
the parameter ao, the terminal bulk speed behaves 
differently. One sees that the terminal speed has 
a maximum for some value of ao (for a given To). 
This is reasonable since for ao ^ and ao — > oo 
we have just one Maxwellian with temperature Tgo 
and for all values of ao between these limits, the 
electron VDF is non-thermal giving rise to a more 
important acceleration because of the velocity fil- 
tration mechanism. In addition, close to these lim- 
its the terminal bulk speed becomes independent 
of To as there is just only one VDF. That makes 
the contour lines to be vertical. 

When using a sum of a Maxwellian core and a 
Kappa halo (for instance with k = 2.5) the con- 
tour plot of the terminal bulk speed is quite simi- 
lar to the previous one as is shown in Fig. 8. The 
main difference is the higher acceleration, which is 
due to the use of the Kappa function. The max- 
ima for a given tq are now displaced to the right 
(to larger values of ao). This is due to the fact 
that for ao — ^ oo the VDF is now just a single 
Kappa accelerating the wind more than a single 
Maxwellian. It is important to note that we can 
obtain very high wind speeds even without using a 
Kappa VDF (Fig. 7). This shows that the acceler- 
ation is not just a consequence of the Kappa func- 
tion, but results from non-thermal distributions, 
as expected. There arc no restrictions on ao and 
To (as in the case of the Kappa VDF where we have 
to take K > 2), but one should constrain these pa- 
rameters by coronal observations or by future in 
situ measurements close to the Sun. 

6. Summary and Final Remarks 

In the present work we have described a col- 
lisionless model of the solar wind acceleration 
assuming non-thermal velocity distribution func- 
tions in the corona. The base altitude of the fast 
wind was taken to be low enough (case of coronal 
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holes), so as to consider the transition from the 
subsonic to the supersonic regime. That needed 
a special resolution method considering a non- 
monotonic potential energy of the protons. An ap- 
proximation regarding the escaping particles rate 
has been used (see Appendix A) in order to obtain 
the electrostatic potential in a self consistent way. 

There are two important results in our work. 
Firstly, we have shown the fundamental role of 
non-thermal electron velocity distributions in ac- 
celerating the wind. The high value of the termi- 
nal bulk speed is not just an artefact of the use 
of Kappa functions. Such speeds can be obtained 
also with a sum of two Maxwellians (a cold and a 
hot one), which is the most commonly used model 
to represent the observed electron distributions. 

Secondly, there is a more important accelera- 
tion of the wind compared to previous exospheric 
models. This is due to the non-monotonic proton 
potential profile that forces some of the protons 
to return back to the Sun and, therefore, reduces 
the escaping proton flux. As a consequence, the 
interplanetary electrostatic potential accelerating 
the wind is enhanced. It is also important to note 
that the terminal bulk speed is anticorrelated with 
the ratio of proton to electron temperatures at the 
base of the wind and therefore there is no need to 
assume large coronal temperatures or additional 
heating of the outer region of the corona in order 
to explain fast wind speeds. 

An inherent property of non-thermal distribu- 
tions is the velocity filtration mechanism (Scuddcr 
1992). This results in an increase of the electron 
temperature within a few solar radii. One problem 
of the present model is the relatively high maxi- 
mum reached by the electron temperature when 
using highly non-thermal electron distributions at 
the cxobasc to produce the fast wind. This prob- 
lem could be due to the approximation used in or- 
der to obtain the potential in a self consistent way. 
A full treatment of the problem with no approx- 
imation would increase the terminal bulk speed, 
but the influence on the electron temperature is 
not known a priori. Another ingredient that might 
influence the results is a possible super-radial ex- 
pansion of the wind. 

An apparent inconsistency of our model is the 
population of trapped electrons. Trapped elec- 
trons can only be produced if collisions decelerate 
ballistic electrons in order to be reflected before 



reaching the cxobasc. Strictly speaking, in the 
present fully coUisionless model, trapped electron 
orbits should not be populated. If so, however, 
the VDF would have a strong discontinuity which 
should be rapidly smoothed out by even a very 
small level of collisions. Hence, we think that it is 
more reasonable - as was done in all previous ex- 
ospheric models - to set this population at a level 
ensuring that the VDF has no discontinuity for 
non-escaping particles. 

The coUisionless nature of our model is an in- 
herent drawback. Recent kinetic simulations of 
the solar wind taking into account binary collisions 
between particles (Landi and Pantellini 2003) sug- 
gest that collisions might be an important in- 
gredient for accelerating the wind to supersonic 
speeds, but this latter work does not consider non- 
thermal electron distributions. The role of colli- 
sions has also been studied by Landi and Pantellini 
(2001) and Dorelli and Scudder (2003) for the so- 
lar corona. The electron heat flux, which plays 
a key role in the solar wind acceleration, seems 
to be essentially determined by the coUisionless 
high-energy tail. This suggests that even though 
the present model neglects collisions, it may cor- 
rectly describe a large part of the physics involved 
in the effect of non- thermal electron distributions. 
In order to better explain the observational prop- 
erties, additional physical ingredients should be 
taken into account, which will be the purpose of a 
future study. 

We are grateful to A. Mangeney, M. Moncu- 
quet and F. Pantellini for stimulating discussions 
and useful comments on the manuscript. We also 
thank the anonymous referee for an amount of 
helpful suggestions that improved this paper. 
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A. Protons accessibility in the E — fx phase space 

In this appendix wc recall some elementary concepts of the resolution method when dealing with protons 
in a non-monotonic potential energy structure. The conservation laws (1) and (2) determine the region where 
the function / is defined as: 

vf>O^E> fj,B{r) + *(r) (Al) 

The relation (Al) defines the line U|| = for each altitude r; the distribution function / is defined only 
above this line, as is shown in Figure 9. From now on we will call this line at the altitude r, as the V limit 
line'. Note that the slope of this line is just the amplitude of the magnetic field B{r) (noted from now on as 
Br). Since in our case, the magnetic field is always decreasing, the sharpest limit line will be the one of the 
reference level tq. Note also that the intercept of an r limit line corresponds to the total proton potential 
energy ^(r) at this altitude (also noted as '^r)- 

Let us now consider the case of a monotonic potential for the protons, i.e. a monotonically decreasing 
potential, which is shown in Figure 10. All limit lines (shown in Fig. 11) for all altitudes r > ro are below 
the ro limit line and they never intersect each other (except for < which is an unphysical case). This is 
due to the fact that both the limit line intercept and its slope are always decreasing. At a given distance 
r, the distribution function / is defined above the corresponding r limit line. But if the particles which are 
present at an altitude r are coming from the exobase ro, their function has also to be defined at the exobase, 
i.e. above the ro limit line. In other words, the only region where the VDF is defined in the E — fi space 
for particles coming from (or getting back to) the exobase is above the ro limit line for all the altitudes. 
The fact that the escaping particle region is defined by only one limit line makes the case of a monotonic 
potential particularly simple. 

The opposite case of a monotonically increasing potential (Fig. 12) is slightly more complicated. The limit 
lines slope is still decreasing, but the intercept is increasing, giving the configuration of Figure 13. Now the 
limit lines for the different altitudes intersect each other. If we consider two altitudes rg and r, with r > ro 
(note that r could be the maximum altitude Vmax in the case of a non-monotonic potential), the VDF / is 
defined, as before, only above the ro limit line and above the r limit line for the altitude r. This means that 
the escaping particles region is not any more defined by only one line, but by two ones that intersect at the 
point (£■*, /u*) with fi* = — (^'r — ^o)/{Br — Bq). This is the a region shown in Figure 13. In the region, 
there are particles which are defined in ro, but are not present in r, that is they are ballistic and can reach 
altitudes up to r and then fall back to the exobase ro. 

In Figures 14 and 15, we can see the same case as before, considering the potential maximum at rmax 
and an intermediate altitude r, with ro < r < rmax- If is the abscissa of the intersection point between 
the ro and r limit lines, we can distinguish two different cases, the first one when /x* > ii^ax^ i-^- the r 
limit line is located above the intersection point {E* ^ji^^^) (Fig. 14) and the opposite one when /x* < fi^ax 
(Fig. 15). In the first case, the escaping region is not defined by only two lines as in Figure 13, because these 
particles have also to be present at altitude r. The consequence is that a new region appears, the region e, 
which contains ballistic particles (between ro and r) and not escaping ones. The totality of ballistic particles 
between ro and r are defined in region 7 + e. The /3 region contains ballistic particles that can reach altitudes 
up to Tmax- We Can therefore see that the escaping particles region a is reduced (by the e region) due to the 
obligatory presence of these particles in the intermediate altitude r between the exobase and Vmax- 

In the second case of an increasing potential (Fig. 15), the escaping region is not modified, but there is 
now a new region {5) with trapped particles between r and rmax, which are defined below the ro limit line 
and therefore do not come from the exobase. It is now evident that definitions of the different species regions 
depend on the potential values at all intermediate altitudes between ro and rmax, as is shown in Figure (16). 
The geometrical definition of these regions is the main numerical difficulty of our problem combined with 
the fact that the potential values are not known but have to be calculated a posteriori using the zero charge 
and current conditions. 

At a given altitude r^, we can consider an elementary region of escaping particles between two consecutive 
points Hi and /X2 and above the r^ limit line (Fig. 17). In order to calculate a moment of the VDF in a given 
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altitude r for the escaping particles, we have to calculate the following integral: 



r rc 

/j(Sa,*a,Ml,M2) = / / 



Fj{E,ii)dE 



djji 



with Ea = + ^'a- Fj IS the appropriate kernel from definitions (3) — (6) with B = Br and = 

For a bi-Maxwellian VDF, these expressions have been calculated by Khazanov ct al. (1998) and will not be 
repeated here. All different integrals Ij have then to be added for all altitudes with a = 0, ...,to. This 
has to be done numerically. 

Up to now, wc have not found a general method for calculating the exact solution of the sclf-consistcnt 
electric potential. Therefore, we have to make an approximation on the escaping particles rate in order 
to find a self-consistent potential. All particles, which can be present at both ro and Vmax, are considered 
as escaping, i.e. wc consider that the region labelled e in Figure 14 corresponds to escaping particles at 
rm , instead of ballistic which do not overcome r. This important approximation is used by Jockers with a 
resulting error in the particle flux claimed to be in general less than 1% (Jockers 1970). However, the errors 
due to this approximation may be much greater with suprathermal electrons. With this approximation 
there are only three unknown parameters, the position of the potential maximum rmax^ the electric potential 
(pEO at the exobase and at the maximum (pEmax- The main problem is that there are only two equations 
(equality of fluxes at the exobase and quasi- neutrality at Vmax), but the problem has a unique solution 
if we suppose the existence of only one maximum, which gives the following constraints: (pEmax < 4'eo-i 

^max > *0=> (pEmax > <pEO - {nipMQG /e){l /vq - l/Vmax) and "^rnax > =^(pEmax > TUpMqG / {ermax) , 

where nip is the proton mass. The technical details about the solution of this system can be found in .Jockers 
(1970) and equally in Lamy et al. (2003). An important point to note is that the approach used by Lamy et 
al. (2003), which consists in integrating the VDF in the velocity space, is leading to exactly the same results 
as our model, using the approximation on the e region. 

In order to verify that the solutions given by the present model are consistent with those of the previous 
models for which the proton potential was monotonic (Maksimovic et al. 1997a), we have calculated the 
electrostatic potential in a range of exobase temperatures including both types of solutions. For simplicity 
we assume for this comparison a single Kappa VDF with Tpo = T^o = To and take ro = 5Rq in order to 
flnd solutions having a monotonic potential in a reasonable temperature range. In Figure 18, we can see the 
electrostatic potential (pEO at the exobase (normalized to kbTo/e) as a function of To for different values of 
K. The full lines correspond to the solutions given by the present model when the proton potential proflle 
is not monotonic. The dashed lines correspond to the solutions obtained with a monotonic proton potential 
proflle (Maksimovic et al. 1997a). These values arc independent of the temperature and represent the lowest 
possible acceleration obtained by exospheric models. The dotted lines belong to the non-monotonic regime as 
well, but the solutions are not known because of a low precision due to the very small difference between (pEo 
and (pEmax- One can see that the solutions obtained in both regimes are mutually consistent as there are no 
discontinuities. The approximation used on the escaping particles rate would have consequences only to the 
very left part of the curves in the Figure 18 for which the acceleration is important and therefore the e region 
would become signiflcant. For the same reasons the approximation made has more important consequences 
for low values of k. In a full treatment of the problem that would not make this approximation, the calculated 
terminal bulk speed would be higher. This is because protons in the e region would be considered as ballistic 
ones (instead of escaping) and would thus need a stronger electrostatic potential in order to escape from the 
gravitational well of the Sun. 
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Fig. 1. — Interplanetary electrostatic potential for 
different values oi k = 2.5, 3.0, 4.0 and 6.0. The 
dashed vertical line indicates Earth's orbit. 
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Fig. 3. Bulk speed profiles for k = 2.5, 3.0, 4.0 
and 6.0. The dashed vertical line indicates the 
Earth's orbit. 
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Fig. 2. — Total proton potential energy for differ- 
ent values of K = 2.5, 3.0, 4.0 and 6.0. The dashed 
vertical line indicates Earth's orbit. The energy is 
normalized to kbTpQ. 
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Fig. 4. — Electron density profiles for k = 2.5, 3.0, 
4.0 and 6.0. The dashed vertical line indicates the 
Earth's orbit. 
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Fig. 5. — Electron temperature profiles for k = 
2.5, 3.0, 4.0 and 6.0. The dashed vertical line in- 
dicates the Earth's orbit. 



18 
1 6 
14 
1 ? 
10 



6 

0.001 



Terminal bulk speed (km/s) 













\ "^""^^^^ 




yy// If r 



































0.010 



0.100 



Fig. 7. — Contours of the terminal bulk speed (at 

lAU) for a sum of two Maxwellians for the elec- 
trons as a function of ao and tq at the exobase. 
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Fig. 6. — Electron temperature profile (full line) 
for a sum of a Maxwellian core and a Kappa halo 
with a = 0.03, r = 5 and k = 2.5. The other 
lines show the contributions of the different parti- 
cle species. The dashed vertical line indicates the 
Earth's orbit. 
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Fig. 8. — As in fig. 7, but for a sum of a Maxwellian 
core and a Kappa halo with k = 2.5. 
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E=|xB(r)+>I'(r) 




Fig. 9. — E — n space for the altitude r. The 
slope of the r limit line is the amplitude of the 
magnetic field -B(r), while its intercept is the total 
proton potential energy ^'(r) at this altitude. The 
distribution function / is defined only above this 
limit line, i.e. in the black region. 




Fig. 10. — Case of an always decreasing proton 
potential energy. 
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Fig. 11. — E — ^ space for the decreasing potential 
of Fig. 10. All r limit lines are found below the rg 
one. The distribution function / is defined in the 
black region for all altitudes r > tq. 




Fig. 12. — Case of an always increasing proton 
potential energy. 




Fig. 13. — E— ^ space for the increasing potential 
of Fig. 12. The r limit lines intersect. Protons that 
can escape from the altitude r are defined in the 
a region. The /3 region corresponds to ballistic 
protons, while the distribution function / is not 
defined in the 7 region. 
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Fig. 14. — E ^ fx space for an increasing proton 
potential energy for altitudes ro < r < rmax, with 
jjL* > jJ'max- The escaping protons region a is re- 
duced by the e one, which contains ballistic pro- 
tons. 




Fig. 16. — E—n space for an increasing proton po- 
tential energy for altitudes ro < r < rmax consid- 
ering the contribution of several altitudes between 

ro and rmax- 



trapped 



Fig. 15. — E — jj, space for an increasing proton 

potential energy for altitudes ro < r < rmax, with 
jjL* < jJimax- The region 5 corresponds to trapped 
particles between r and Vmax- 




Fig. 17. — An example of an elementary region of 
escaping particles between two consecutive points 
/ii and /U2 and above the limit line. 
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To (lO^K) 



Fig. 18. — Electrostatic potential in function of 
the exobase temperature for different values of k. 
Full lines represent solutions given by the present 
model (non-monotonic proton potential). Dashed 
lines correspond to previous exosphcric models 
with a monotonia proton potential profile. Dotted 
lines show the range of a low numerical precision 
as explained in the text. 



